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29 Abstract 

30 A 3-D cloud scale chemical transport model that includes a parameterized source 

3 1 of lightning NO x based on observed flash rates has been used to simulate six midlatitude 

32 and subtropical thunderstorms observed during four field projects. Production per 

33 intracloud (Pic) and cloud-to-ground (Pcg) flash is estimated by assuming various values 

34 of Pic and Pcg for each storm and determining which production scenario yields NO x 

35 mixing ratios that compare most favorably with in-cloud aircraft observations. We obtain 

36 a mean Pcg value of 500 moles NO (7 kg N) per flash. The results of this analysis also 

37 suggest that on average, Pic may be nearly equal to Pcg, which is contrary to the common 

38 assumption that intracloud flashes are significantly less productive of NO than are cloud- 

39 to-ground flashes. This study also presents vertical profiles of the mass of lightning NO x 

40 after convection based on 3-D cloud-scale model simulations. The results suggest that 

41 following convection, a large percentage of lightning NO x remains in the middle and 

42 upper troposphere where it originated, while only a small percentage is found near the 

43 surface. The results of this work differ from profiles calculated from 2-D cloud-scale 

44 model simulations with a simpler lightning parameterization that were peaked near the 

45 surface and in the upper troposphere (referred to as a “C-shaped” profile). The new 

46 model results (a backward C-shaped profile) suggest that chemical transport models that 

47 assume a C-shaped vertical profile of lightning NO x mass may place too much mass near 

48 the surface and too little in the middle troposphere. 

49 
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49 1 . Introduction 

50 The oxides of nitrogen, (NO + NO2 = NO x ), are important O3 precursors in the 

5 1 troposphere. Of the major sources of NO x in the troposphere, lightning remains the 

52 source with the greatest uncertainty and is particularly important because it produces NO x 

53 in the middle and upper troposphere where NO x is longer lived and can be more efficient 

54 at producing ozone than in the boundary layer. The representation of lightning NO x 

55 (LNO x ) in three-dimensional (3-D) regional and global chemical transport models 

56 (CTMs) is critical to the models’ representation of ozone and other species such as OH 

57 [e.g. Stockwell et al., 1999, Labrador et al., 2004], Labrador et al. [2005] found that both 

58 the magnitude of the global LNO x source strength and its vertical distribution can 

59 substantially affect tropospheric trace gas concentrations in a global CTM. IPCC [2007] 

60 suggests a global LNO x source of 1 .1-6.4 Tg N yr" 1 based on the work of Boersma et al. 

61 [2005]. However, Schumann and Huntrieser et al. [2007] have comprehensively 

62 summarized estimates of global LNO x production and found the best estimate to be 5±3 

63 Tg N yr’ 1 . 

64 In order to adequately represent the LNO x source in global or regional models, the 

65 geographic distribution of flashes, average production of NO per flash, and vertical 

66 distribution of LNO x following convection must be specified. A variety of schemes have 

67 been used to specify the horizontal distribution of flashes (e.g. using variables such as 

68 cloud-top height, upward cloud mass flux, convective preciptation, and CAPE as 

69 predictors). Allen and Pickering [2002] have evaluated their use in global 3-D CTMs. 

70 The production of NO per flash has also been examined using laboratory' experiments. 

71 theoretical assumptions regarding the physics of lightning flashes, and observations 

72 obtained during field projects. Despite these efforts, a great deal of uncertainty remains 


3 



73 regarding NO production on a per flash basis, as well as the relative production by 

74 intracloud (IC) and cloud-to-ground (CG) flashes. On the basis of previous studies, 

75 which suggested that IC flashes were less energetic than CG flashes (e.g. Holmes et at., 

76 1971), many studies of lightning NO x production have assumed that Pic is less than Pec,. 

77 Price et al. [1997] assumed that a CG flash produces approximately 1100 moles of NO 

78 and that an IC flash was one tenth as productive of NO as a CG flash in estimating global 

79 lightning NO* production. 

80 However, a number of more recent studies have suggested that Pic may be nearly 

81 as great as Pcg- Gallardo and Cooray [1996] suggested that IC flashes may dissipate 

82 nearly as much energy as CG flashes and therefore Pic may be on the order of Pcg- 

83 Supporting the Gallardo and Cooray [1996] hypothesis, a two-dimensional (2-D) cloud- 

84 scale modeling study by DeCaria et al. [2000] suggested that the Pic/Pcg ratio is likely 

85 between 0.5 and 1 .0 and a 3-D modeling analysis by DeCaria et al. [2005] narrowed this 

86 range to between 0.75 and 1.0. A study by Fehr et al. [2004] used a 3-D cloud-scale 

87 model simulation of a storm observed over Germany, and by comparing with lightning 

88 and aircraft observations, concluded that on average, an IC flash produced 40% more NO 

89 than a CG flash. Zhang et al. [2003] argued that IC flashes may dissipate 50-100% as 

90 much energy as CG flashes because they have a large number of return strokes which 

91 where neglected in the Price et al. [1997) calculations. Rahman et al. [2007] found that 

92 results from experiments with rocket-triggered lightning suggested that NO x production 

93 was associated with relatively long duration continuing currents, which may be greater in 

94 IC flashes than in CG flashes, leading to the implication that IC flashes may produce as 

95 much or possibly more NO x per flash than CG flashes. Recent estimates of energy 

96 dissipated by IC and CG flashes using electrical potential and charge density 
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97 measurements by Maggio et al. [2008] also suggest IC flashes may dissipate as much or 

98 more energy than CG flashes. 

99 Pickering et al. [1998] presented vertical profiles of lightning NO x for use in 3 D 

1 00 CTMs based on the results of 2-D cloud-resolving model simulations of seven convective 

101 events. These simulations assumed the production scheme of Price et al. [1997]. NO, 

1 02 produced by CG flashes was distributed in the simulated storms from the surface to the - 

103 1 5°C isotherm while NO x produced by IC flashes was distributed from the -1 5°C 

1 04 isotherm to the cloud top. Average profiles of LNO x mass computed for the midlatitude 

1 05 continental, tropical continental, and tropical marine regimes showed peaks in mass near 

1 06 the surface and in the upper troposphere, leading many CTMs to adopt a C-shaped 

1 07 vertical distribution of LNO x mass. 

1 08 This paper has two primary objectives: (1) to summarize the results of 3-D cloud- 

1 09 resolved storm simulations yielding estimates of NO production per flash and (2) to 

110 update the vertical LNO x profile information of Pickering et al. [1998]. Six storms from 

1 1 1 field projects conducted in Germany, Colorado, south Florida, and Kansas/Oklahoma 

1 12 have been simulated using a 3-D cloud-scale chemical transport model (CSCTM) that 

1 1 3 includes a parameterized source of LNO x . LNO x production per flash is estimated 

1 1 4 individually for five of the six storms. In this study we present vertical distributions of 

1 1 5 LNO x calculated from each of these 3-D simulations as well as average vertical profiles 

116 for the midlatitude continental and subtropical events, which can be applied in regional 

1 1 7 and global CTMs. Section 2 describes the methodology used in these studies, while 

1 1 8 section 3 presents results from the individual storm simulations. Section 4 discusses the 

1 1 9 results of the storm case studies and their application to global models and remote 

1 20 sensing. Section 5 presents conclusions that may be drawn from this work. 
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121 2. Methodology 

122 Storms from the Cirrus Regional Study of Tropical Anvils and Cirrus Layers - 

1 23 Florida Area Cirrus Experiment (CRYSTAL-FACE; Ridley et al., 2004; Lopez et al., 

124 2006), European Lightning Nitrogen Oxides Project (EULINOX; Huntrieser et al, 2002). 

125 Stratosphere Troposphere Experiment: Radiation, Aerosols and Ozone (STERAO; Dye et 

126 al., 2000), and Preliminary Regional Experiment for STORM (PRE-STORM; Rutledge 

127 and MacGorman, 1988) field projects were simulated. With the exception of the PRE- 

128 STORM event, all of these storms featured measurements of chemical and 

1 29 meteorological properties by research aircraft at anvil levels. The time, location, and 

130 peak current of CG lightning occurrences in all storms were recorded by ground-based 

1 3 1 systems, and during STERAO and EULINOX, total lightning activity (IC + CG) was 

1 32 mapped by a VHF interferometer. In addition, all experiments included extensive 

133 satellite and radar observations of storm development and evolution. 

134 The dynamical evolution of each storm was simulated using a cloud-resolving 

135 model and the temperature, wind, and hydrometeor fields were then used to drive the 

136 offline CSCTM. For each storm, various LNO x production per flash scenarios were 

137 simulated and model results were compared with in-cloud aircraft observations of NO x 

138 through use of mean profiles, column NO x mass, and probability distribution functions of 

139 simulated and observed NO x to determine the most appropriate scenario. At the end of 

140 the CSCTM simulation, the mass of N fixed by lightning was calculated at each model 

141 level, and the percentage of the total mass of LNO x was calculated for 1-km layers. 

142 A detailed description of the CSCTM is found in DeCaria et al. [2005], In this 

143 version of the model, LNO x production is computed using observed IC and CG flash 

1 44 rates and a specified scenario of Pic and Pco to calculate the mass of NO injected into the 
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145 cloud per time step. Both laboratory experiments (Wang et al., 1998) and theoretical 

146 considerations (Price et al., 1997) have pointed to a strong dependence of LNO x 

147 production on stroke peak current. Therefore, we have made initial estimates of Pco 

148 using these relationships because peak current data for return strokes are available from 

149 ground-based network observations of CG flashes. However, this assumption contributes 

150 to uncertainty in the results because recent work with rocket-triggered lightning (Rahman 

151 et al., 2007) has suggested that return strokes are not the primary NO-producing phase of 

1 52 a lightning flash. Field experimental results from TROCCINOX (Tropical Convection, 

1 53 Cirrus and Nitrogen Oxides Experiment) in Brazil (Huntrieser et al., 2008) have 

1 54 suggested a more minor dependence on stroke peak current. 

1 55 The NO produced by CG flashes is distributed unimodally in the vertical, while 

1 56 the NO produced by IC flashes is distributed bimodally based on the vertical distributions 

157 of very high frequency (VHF) sources of IC and CG flashes presented in MacGorman 

158 and Rust [1998]. The two modes correspond to the two main charge centers in a typical 

159 thunderstorm cloud. The vertical distributions of IC and CG lightning channels used in 

160 the model are shown in Figure 2 of DeCaria et al. [2005]. The modes of the IC 

161 distribution are nominally set at the heights of the -15°C and -45°C isotherms. However, 

162 the vertical distribution of IC channels is modified as necessary in some storms by 

1 63 changing the upper mode isotherm to account for a higher or lower cloud top and to 

164 match the upper tropospheric peak in aircraft observations of NO x . At each model level. 

165 the lightning NO is distributed uniformly to all grid cells within the 20 dBZ contour 

1 66 computed from simulated hydrometeor fields. A passive version of the CSCTM includes 

167 only the transport of tracer species and production of lightning NO x . In this version, 

168 three types of NO x (pre-existing NO x , CG LNO x , and IC LNO x ) are transported without 
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169 any chemical reactions and the results are used to estimate average NO x produced by 1C 

170 and CG flashes. Various production scenarios are specified and the model results 

171 compared with in-cloud aircraft observations to determine a first estimate of the most 

172 appropriate LNO x production scenario for the storm in question. 

1 73 A version of the model that includes O 3 photochemistry is used to separate NO x 

174 into NO and NO 2 and obtain a final and best estimate of NO production by IC and CG 

1 75 flashes. The chemical mechanism for the STERAO storms is described in DeCaria et al. 

176 [2005]. In the CRYSTAL-FACE and EULINOX simulations, the chemical scheme is 

1 77 identical to that described in DeCaria et al. [2005] with the exception that in these cases, 

178 isoprene and propene chemistry was included. Initial profiles of ozone were constructed 

179 using out-of-cloud aircraft observations, and when necessary to fill gaps, climatological 

1 80 ozone profiles appropriate for the latitude of the storm in question. The CSCTM include 

181 a simple scheme to represent the influence of clouds on photolysis rates based on the 

1 82 work of Madronich [ 1 987] and fully described in DeCaria et al. [2005], Clear sky 

183 photolysis rates are calculated following Stamnes et al. [1988] using observed column 

1 84 ozone amounts measured by the Total Ozone Mapping Spectrometer (TOMS) or, in the 

1 85 case of the STERAO storms, ozonesondes launched from Boulder, Colorado. 

1 86 Perturbations to clear sky photolysis rates are determined by cloud thickness with very 

1 87 thick clouds (deeper than 5 km) resulting in photolysis rates multiplied by a factor of 2 

1 88 above the cloud and 0. 1 below the cloud base with a linearly interpolated degree of 

1 89 enhancement at intermediate altitudes. 

190 3. Results 

191 3.1 Subtropical Events 
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192 The CRYSTAL-FACE experiment was conducted over South Florida in July, 

1 93 2002 and two storms investigated during the campaign were simulated. The July 29 

194 CRYSTAL-FACE storm was simulated using the NASA Goddard version of the non- 

195 hydrostatic PSU/NCAR (MM5) mesoscale model (Tao et al, 2003b) with a horizontal 

196 resolution of 2 km and vertical resolution of 0.5 km. At 1700 UTC (1300 LT) on July 29 , 

1 97 2002, a powerful thunderstorm developed along the west coast of Florida near Fort 

1 98 Myers. While the storm intensified and moved north along the coast, the area in and 

199 above the anvil was sampled by the NASA WB-57 research aircraft from 1845 to 201 1 

200 UTC at altitudes ranging from 12.5 to 13.8 km. The coastal convection later merged with 

201 convection originating near Lake Okeechobee. Figure 1 shows the time series of CG 

202 flashes recorded by the National Lightning Detection Network (NLDN) from 1700 to 

203 2300 UTC. The July 29 storm was an exceptionally strong lightning producer with 4168 

204 CG flashes recorded during this period. Maxi m um CG flash rates exceed 30 flashes per 

205 minute. These lightning flash rates and their relationship to NO observations from the 

206 WB-57 aircraft have been discussed by Ridley et al. [2004], 

207 Simulated convection along the coast began earlier than observed by 

208 approximately three hours, but a number of storm features were successfully reproduced, 

209 including cloud top height and the direction of storm movement. The transport of CO 

210 (initial boundary layer maximum) and O 3 (initial stratospheric maximum) were calculated 

211 by the CSCTM and compared with in-cloud observations using probability distribution 

212 functions (pdfs). The model overestimated the frequency of values in the upper end of 

213 the frequency distributions for both CO and O3 at anvil levels, suggesting that both 

214 upward and downward transport in the model may be too strong. However, the intense 

2 1 5 lightning activity in the thunderstorm and extremely elevated NO x mixing ratios (up to 1 0 
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216 ppbv) observed in the anvil indicate that the lightning NO x source in this storm was 

217 extremely strong. As a result, errors in transport are unlikely to significantly affect the 

218 estimate of LNO x production per flash in this case. 

2 1 9 Because the NLDN only recorded the occurrence of CG flashes during 

220 CRYSTAL-FACE, IC flashrates were estimated for the July 29 storm. NLDN flashes 

221 with positive peak current less than 10 kA are thought to be IC flashes [Cummins et al., 

222 1 998], These flashes were removed from the NLDN observations to estimate CG flash 

223 rates. The percentage of recorded flashes with positive peak current less than 10 kA was 

224 calculated during both the July 29 storm and the month of July as a whole. The 

225 percentage of such flashes in the July 29 storm was larger by a factor of 2.5 than the 

226 percentage for the month of July which may indicate an enhancement of the IC/CG ratio 

227 in the July 29 storm over the climatological value. To account for this enhancement, the 

228 south Florida climatological value for the IC/CG ratio of 2 from Boccippio et al. [2001] 

229 was multiplied by 2.5 to estimate an IC/CG ratio of 5. The upper mode of the vertical 

230 distribution of IC flash channel segments was set to -45°C. 

23 1 P cq was estimated to be approximately 590 moles NO based on observed mean 

232 peak current (19 kA) and a relationship between peak current and energy dissipated from 

233 Price et al. [1997]. Various values of the Pic/Pcg ratio were simulated and the results 

234 compared with aircraft NO x observations (observed NO + estimated NO 2 ). NO 2 was 

235 estimated using NO and O 3 observations and the photostationary state assumption 

236 because no direct measurements of NOi were available. Following the work of 

237 Madronich ( 1 987), who demonstrated that actinic flux and photodissocation could 

238 increase within clouds, NO 2 mixing ratios were calculated assuming both clear sky 

239 photolysis rates and rates enhanced by a factor of 2. The observed column mass of N in 


10 



240 NO x was calculated for the 1 km thick layer extending from 12.25 to 13.25 km by 

241 computing layer mean in-cloud NO x mixing ratios. Assuming that photolysis rates were 

242 enhanced in the cloud yielded a column mass of 6.6 x 10' 4 g N m' 2 while assuming that 

243 photolysis rates were unaltered by the cloud yielded a column mass of 7.0 x 10' 4 g N nf 2 

244 Table 1 shows the sensitivity of CSCTM-calculated column mass to the assumed Pic/Pec 

245 ratio. Assuming that an IC flash is on average 60% as productive of NO as a CG flash 

246 yields the most favorable comparison with column mass of nitrogen estimated from 

247 aircraft observations when an IC/CG ratio of 5 is assumed. 

248 Figure 2 shows the pdfs of observed and simulated in-cloud NO x assuming that 

249 Pic equals 50 and 60% of Pcg and an IC/CG ratio of 5. When compared with observed 

250 pdfs at both 12.5 and 13 km, the production scenario that assumes a Pic/Pcg ratio of 0.5 

25 1 compares more favorably than assuming a Pic/Pcg ratio of 0.6. The pdfs show that the 

252 underestimation of the simulated column mass assuming a Pic/Pcg ratio of 0.5 results 

253 from the inclusion of gridcells with background NO x concentrations as well as the 

254 underestimation of the distribution at 13 km between 2 and 6 ppbv. Based on the 

255 comparison of simulated and observed pdfs and computed column mass, it is estimated 

256 that in the July 29 CRYSTAL-FACE storm an IC flash on average produced 50-60% as 

257 much NO as a CG flash while an average CG flash produced 590 moles NO. The 

258 assumption that an IC flash is one tenth as productive of NO as a CG flash significantly 

259 underestimated NO x when compared to column mass estimates calculated from in-cloud 

260 aircraft observations. When chemical reactions were simulated in addition to lightning 

261 NO x production and convective transport, a loss of NO x occurred due to conversion to 

262 reservoir species. The simulated column mass of N in NO x decreased by approximately 

263 5% from the values shown in Table 1 . The Pic/Pcg ratio of 0.6 was selected for the 
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264 subsequent calculations because it continued to provide the best comparison with the 

265 column mass estimates derived from aircraft. It should be noted that the production 

266 scenario discussed above was deduced assuming that the IC/CG ratio in this particular 

267 storm was greater than the climatological IC/CG ratio for south Florida by a factor of 2.5. 

268 Because many more weak positive flashes (which are believed to be IC flashes) were 

269 recorded in this storm than was typical for the South Florida area during the month of 

270 July, it is likely that the IC/CG ratio was elevated above the climatological value. If 

271 instead the climatological IC to CG ratio of 2 is assumed, an IC flash must produce 50% 

272 more NO than a CG flash to match aircraft observations. 

273 It should be noted that the assumptions of IC/CG ratio, Pcg value, and Pic/Pco 

274 ratio contribute to uncertainty in the production scenario estimate. Figure 3 shows all the 

275 possible production scenarios that yield NO x mixing ratios which would match observed 

276 column mass when IC/CG ratios of 2 or 5 are assumed. The symbols on the curves show 

277 our best estimates (Pcg = 590 moles/flash with Pic/Pcg at 0.6 and 1.5) based on the 

278 observed mean peak current. Most literature estimates of Pcg range from 200 to 1 1 00 

279 moles NO per flash. For this range of Pcg values, the Pic/Pcg ratio could range from 0.3 

280 to 4 in the July 29 CRYSTAL -FACE storm. If Pic is assumed to equal Pcg, then Pcg 

28 1 must be between 350 and 850 moles NO per flash depending on the number of IC flashes 

282 which are assumed to have occurred. This analysis puts bounds on the uncertainty 

283 associated with assuming a single production scenario. If the IC/CG ratio was known to 

284 be 5 for this storm, the uncertainty would be reduced. Pic/Pcg would be in the 0.3 to 1 .8 

285 range (assuming Pcg is between 200 and 1 100 moles NO per flash). The assumption that 

286 IC and CG lighting are equally productive of NO on a per flash basis would yield a 

287 production of 350 moles NO per flash. This analysis demonstrates the need for 
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288 observations of IC flashes in future experiments and more definitive information on the 

289 relationship between peak current or other electrophysical variables and LNO x 

290 production. 

291 A vertical cross section of simulated NO x through the core of the coastal storm at 

292 the end of the 240 minute simulation is shown in Figure 4 (assuming Pcg = 590 moles 

293 NO, Pic = 354 moles NO, and an IC/CG ratio of 5). Maximum NO x mixing ratios exceed 

294 1 1 ppbv in the convective plume extending west from the Florida coast in the 10 to 12 

295 km region. Figure 5 shows the vertical profile of the lightning NO x mass as it is 

296 introduced into the model domain, as well as its profile following convection from the 

297 same simulation as Figure 4. The bimodal distribution of the profile of injected LNO x 

298 mass reflects the bimodal distribution of IC flash channel segments calculated by the 

299 CSCTM. The upper mode peak is initially smaller in magnitude because the LNO x is 

300 distributed with a dependence on pressure as described in DeCaria et al. [2005], Upward 

301 transport during the storm results in LNO x originally introduced into the model at 

302 altitudes near the lower mode peak residing near the top of the cloud following 

303 convection. 

304 On July 16, 2002, an isolated convective system developed northwest of Miami 

305 shortly after 1900 UTC (1500 LT) and was investigated as part of the CRYSTAL-FACE 

306 project. Over the next few hours, the storm moved west across the Florida peninsula arid 

307 was extensively sampled by the WB-57 from 1936 to 2306 UTC at altitudes ranging from 

308 9 to 15.5 km AGL. The storm was simulated using the Advanced Regional Prediction 

309 System (ARPS) described in Xue et al. [2000; 2001] with a horizontal resolution of 2 km 

3 1 0 and vertical resolution varying from 25 m near the surface to 0.5 km near the top of the 

3 1 1 model domain at approximately 25 km. A number of different types of observations. 
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312 including radar reflectivity, were assimilated into the simulation. The simulated temporal 

3 1 3 evolution of the storm matched observations well, as did storm size. The CSCTM was 

314 used to calculate the transport of tracer species CO and O 3 and a comparison of model 

3 1 5 results with anvil aircraft observations showed that the model adequately represented 

3 1 6 transport within the storm. 

3 1 7 Figure 6 shows the time series of flash rates recorded by the ground-based NI.DN 

3 1 8 from 1 900 to 2300 UTC. In contrast to the July 29 CRYSTAL-FACE storm, the July 1 6 

319 storm was a relatively weak lightning producer with only 301 CG flashes recorded during 

320 this period and a maximum CG flash rate of 9 flashes per minute. The percentage of 

321 weak positive flashes (peak current < lOkA) was calculated during the July 16 storm and 

322 was over 20 times greater during the storm than for the month of July as a whole. 

323 Because the flash rate during the storm represented too small a sample to reliably make a 

324 large adjustment, IC flash rates were estimated from the observed CG flash rates and the 

325 climatological IC to CG ratio of 2 for south Florida from Boccippio et al. [2001], Based 

326 on the mean peak current of CG flashes recorded by the NLDN for this storm (23 kA), 

327 Pcg was estimated to be 700 moles of NO. In this simulation, the upper mode of the IC 

328 flash channel distribution was set to -60°C. Several different values of the Prc/Pco ratio 

329 were simulated and the results compared with observations. The assumption that on 

330 average, an IC flash produces 75% as much NO as a CG flash yielded the most favorable 

33 1 comparison with the column mass of N in NO x estimated from observations when 

332 chemical reactions were not simulated. The simulation of chemistry led to a decrease in 

333 NO x mixing ratios due to conversion of NO x to reservoir species such as UNO 3 and PAN. 

334 In order to match aircraft observations, a Pic/Pcg ratio of 0.9 was needed. Figure 7 shows 

335 a vertical cross-section of NO x calculated assuming this production scenario taken 
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336 through the core of the storm after 180 minutes of simulation. Maximum NO x mixing 

337 ratios exceeding 5 ppbv are found at 6.5 km while lower concentrations are found in the 

338 upper part of the core and anvil. 

339 Figure 8 shows the percentage of the mass of N in LNO x in each kilometer-deep 

340 layer after both the July 16 and 29 CRYSTAL-FACE storms. Following the July 29 

341 storm, the maximum in the vertical mass distribution is found ait anvil levels ( — 10 — 11 

342 km). In the case of the July 16 storm, the maximum is found in the 6-7 km layer, 

343 coincident with the lower mode of the vertical distributions of the LNO x source in the 

344 CSCTM. A smaller peak is found in 9-10 km layer, the height of the upper mode of the 

345 distribution, and another peak at 12-13 km, near the top of the cloud. These vertical 

346 LNO x mass distributions are based on the assumed IC to CG ratios of 2 and 5 for the July 

347 1 6 and 29 storms, respectively. Because IC to CG ratios may be highly variable and were 

348 estimated for these simulations, the results of a sensitivity test of the assumption of IC to 

349 CG ratios of 2 and 5 are shown for the July 16 storm in Figure 9a and for the July 29 

350 storm in Figure 9b. For both storms, assuming an IC to CG ratio of 2 results in a slightly 

351 larger percentage of LNO x mass residing in the lower and middle portion of the cloud and 

352 a smaller percentage of LNO x mass near the cloud top (12-14 km for the July 16 case and 

353 10-12 km for the July 29 storm). 

354 3.2. Midlatitude Continental Events 

355 The STERAO field project was conducted in June and July, 1996 over 

356 northeastern Colorado to study lightning NO x production and convective transport of 

357 chemical species. CG lightning activity was monitored by the NLDN, while total 

358 lightning activity, including both CG and IC flashes, was mapped by a VHF 

359 interferometer. Two storms observed during the STERAO field project were simulated 
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360 and lightning NO x production per flash estimated in each. A similar analysis was 

361 conducted for one storm from the EULINOX field project conducted in southern 

362 Germany in July 1998. The fourth midlatitude storm simulated was from the 1985 

363 PRESTORM experiment conducted over Kansas/Oklahoma. 

364 All four midlatitude storms were simulated using the 3-D Goddard Cumulus 

365 Ensemble (GCE) model, which is described in Tao and Simpson [1993] and Tao et al. 

366 [2003a], A description of the GCE simulation of the July 12 STERAO storm is contained 

367 in Stenchikov et al. [2005], while the results of the chemical transport model simulation 

368 are found in DeCaria et al. [2005]. Between 2000 and 2400 UTC (1400 and 1800 LT), 

369 188 CG flashes and 2121 IC flashes were recorded by the NLDN and interferometer, 

370 respectively. Flash counts for the STERAO storms only include interferometer flashes 

371 with duration greater than 100 ms. Whether or not short duration flashes, which would 

372 not be recorded by some lightning detection systems, are productive of NO is an open 

373 question. DeCaria et al. [2005] estimated that Pcg was 460 moles NO based on the mean 

374 peak current of 15 kA. Assuming that an IC flash was equally as efficient at producing 

375 NO as a CG flash yielded the best comparison with observations in terms of mean 

376 vertical profile shapes while a Pic/Pcg ratio between 0.75 and 1 compared most favorably 

377 with observed column mass. 

378 On July 10, 1996, a multi-cellular thunderstorm organized in a northwest to 

379 southeast line developed near the Wyoming-Nebraska border at approximately 2100 UTC 

380 (1500 LT) and was observed as part of the STERAO campaign [Dye et al., 2000; 

381 Skamarock et al., 2000; 2003]. Research aircraft investigated the storm from 2237 to 

382 0105 UTC. After 0115 UTC, the storm became unicellular with supercell chracteristics 

383 IC lightning activity dominated throughout the lifetime of the storm with only 77 CG 
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384 flashes and 3223 IC flashes recorded between 2201 and 0206 UTC. The University of 

385 North Dakota Citation aircraft observed maximum NO mixing ratios of 1 ppbv 

386 approximately 60 km downstream of the storm cores [Dye et al., 2000]. 

387 The GCE model with a horizontal resolution of 2 km and 0.5 km vertical 

388 resolution was used to simulate the July 10 STERAO storm. The GCE and CSCTM 

389 simulations were included in a cloud-chemistry model intercomparison study described in 

390 Barth et al. [2007], Both the magnitude and height of the simulated peak updraft velocity 

391 compared well with observations. Convective transport was evaluated by comparing 

392 simulated CO and O 3 mixing ratios with observations obtained during two across-anvil 

393 transects. While the simulation underrepresented transport during the first transect, 

394 mixing ratios in the later transect were reproduced well by the CSCTM simulation, 

395 indicating that convective transport was adequate after an initial spin-up period. 

396 Simulated storm speed and cloud top heights compared favorably with radar 

397 observations. For the simulations of the July 10 STERAO and July 21 EULINOX 

398 storms, the parameterization of LNO x production was modified to simulate individual 

399 flashes. The method is fully illustrated in Ott et al. [2007]. Each flash is simulated by 

400 selecting grid cells at random from an area of the domain centered just downwind of the 

401 maximum updraft, which is similar in size to the area where lightning flashes typically 

402 occurred. The downwind location was chosen based on plots of radar reflectivity and 

403 flash location in Holler et al. [2000] and Dye et al. [2000] that show the majority of 

404 lightning activity occurring slightly downwind of the core updraft region. The number of 

405 grid cells included in a flash at each level is determined by vertical distributions adapted 

406 from the Gaussian distributions calculated in the original bulk method of LNO x 

407 parameterization described in DeCaria et al. [2005]. Simulations of both the July 10 
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408 STERAO and July 21 EULINOX storms using both the original and modified LNO x 

409 parameterizations show that the difference in parameterization produces little change in 

410 the vertical distribution of LNO x following convection. The upper mode of the IC flash 

41 1 channel distribution was set to -50°C for the July 10 STERAO storm. 

412 In the July 10 STERAO storm, Pcg was estimated to be approximately 390 moles 

413 based on the observed peak current (13 kA) of CG flashes. Assuming a Pic/Pcg ratio of 

414 0.5 matched in-cloud aircraft observations well when chemical reactions were not 

415 considered. When chemical reactions were simulated, the Pic/Pcg ratio was increased to 

416 0.6 to match observations. Figure 10 shows a vertical cross section of calculated NO x 

4 1 7 through the core of the southernmost cell in the line of thunderstorms at the end of the 

418 180 minute simulation. NO x mixing ratios exceed 2.8 ppbv at 3.5 km above the surface 

419 but are only ~1 ppbv in the anvil region. 

420 The EULINOX field project was conducted in central Europe in June and Jub 

42 1 1 998 in order to study lightning NO x production over Europe. In addition to two research 

422 aircraft, the project included observations of CG lightning activity from the ground-based 

423 lightning detection system known as BLIDS (Blitz Informationsdienst von Siemens). 

424 One storm from the EULINOX project, which occurred on July 21, 1998 over southern 

425 Germany, was simulated using the GCE model. The storm developed as a single cell and 

426 after an initial period of intensification split into two distinct cells. The northernmost cell 

427 became multicellular in structure and was observed to decay rapidly, while the southern 

428 cell strengthened and developed supercell characteristics. The southern cell produced 

429 360 CG and 2565 IC flashes between 1640 and 1900 UTC (1840 and 2100 LT), while the 

430 northern storm produced 289 CG and 815 IC flashes. The GCE simulation succeeded in 

43 1 reproducing a number observed storm features, including the cell-splitting. A full 
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432 description of the GCE and CSCTM simulations is provided in Ott et al. [2007]. Based 

433 on the mean peak current (12 kA) of CG flashes, Pcg was estimated to be 360 moles NO. 

434 The upper mode of the IC flash channel distribution was set to -45°C in this simulation. 

435 Several different values of Pic were simulated, and the scenario in which Pic is equal 10 

436 Pcg was found to compare most favorably with observations, though observed column 

437 mass was underestimated by 10%. In order to match observed column mass, the Pic/P< c; 

438 ratio needed to be increased to 1.15. 

439 In addition to the STERAO and EULINOX storms, the June 10-11 squall line 

440 observed during the PRE-STORM project was simulated using the GCE model. The 

441 June 10-11 squall line has been documented extensively [e.g. Johnson and Hamilton, 

442 1988; Rutledge et ah, 1988] and has previously been simulated using the 2-D version of 

443 the GCE [Tao et al., 1993]. In this case, the horizontal resolution of the 3-D GCE was 

444 1.5 km and vertical resolution varied from approximately 0.25 km near the surface to 

445 slightly more than 1 km near the top of the domain at 2 1 .4 km. During PRE-STORM, the 

446 occurrence of CG lightning flashes was recorded by the National Severe Storms 

447 Laboratory's Lightning Location Network. Based on a time series of positive and 

448 negative CG flash rates from Nielsen et al. [1994], approximately 6500 CG flashes 

449 occurred during the storm's lifetime. Observations of total lightning activity were 

450 unavailable, so the climatological IC/CG ratio of 3 for the region [Boccippio et al, 200 ' | 

45 1 was assumed to estimate IC flash rates. The upper mode of the IC flash channel 

452 distribution was set to -45°C. No observations of the chemical environment of the squall 

453 line anvil are available. Therefore, it was not possible to estimate a production scenario 

454 for IC and CG flashes as in the other five storms. Instead, the average value of Pcg 
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455 calculated over the five other storms (~500 moles NO per flash) was used along with an 

456 estimate of Pic that was 85% of Pcg or 425 moles of NO. 

457 Figure 1 1 shows the vertical distribution of the mass of N in LNO x for the four 

458 midlatitude continental storms. The distributions for the four storms all reflect the double 

459 peaked distribution of LNO x produced by IC flashes in the model. There is variation 

460 between the simulations in the dominant mode of the lightning distributions. In the 

461 EULINOX and PRESTORM storms whose IC to CG ratios were on average 5 and 3, 

462 respectively, a higher percentage of LNO x mass is found near the height of the lower 

463 mode of the IC distribution, which is also the mode of the CG distribution. In the July i 0 

464 and July 12 STERAO storms, which had average IC to CG ratios of 33 and 8 during the 

465 time periods simulated, a greater percentage of LNO x mass resides near the height of the 

466 upper mode of the IC vertical distribution following convection. In addition to the IC to 

467 CG ratio, the dominance of the modes is likely affected by storm dynamics and the 

468 timing of IC and CG flashes in relation to the evolution of the storm. 

469 4. Applications for large-scale models and remote sensing 

470 4.1 Vertical distribution of lightning NO x mass 

471 We have presented results from the simulations of six thunderstorms using the 3 - 

472 D CSCTM. Figure 12a shows the average vertical distribution of the mass of N in LNO x 

473 calculated by averaging the case studies in the subtropical regime and Figure 12b shows 

474 the average vertical distribution for storms in the midlatitude regime overlaid with the 

475 vertical distribution calculated by Pickering et al. [1998]. Both plots are overlaid with 

476 smooth curves fit to the regime average. Table 2 lists the percentages of LNO x mass in 

477 each 1 -km layer taken from the smoothed curves. In both regimes on average, only a 

478 small percentage of LNO x resides in the boundary layer following the convective event. 
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479 A greater percentage of LNO x remains in the middle and upper troposphere where the 

480 LNO x was originally produced. These “backward C-shaped” average vertical 

481 distributions are in marked contrast to the C-shaped profiles presented in Pickering et al. 

482 [1998] based on 2-D cloud-resolving model simulations where a significant percentage oi 

483 LNO x mass was transported to the boundary layer and relatively little LNO x mass was 

484 found between 1 .5 and 6.5 km after convection concluded. Our results are similar to 

485 those from recent 3-D cloud models with explicit electrophysics. For example, Zhang et 

486 al. [2003a] found a maximum in NO x at mid-levels along with a secondary maximum at 

487 anvil levels indicating that a large portion of the LNO x remains within the cloud near the 

488 levels of its production. It should also be noted that our simulations indicate a small 

489 percentage of mass (~1%) remains above the tropopause following convection, which 

490 results from overshooting transport of lightning NO x produced at lower altitudes. 

491 Individual modeling groups who implement these profiles may want to consider scaling 

492 these profiles so that mass is not directly injected into the stratosphere because most 

493 convective parameterizations used in larger-scale models do not produce overshooting 

494 cloud tops. Implementation of the profiles should involve scaling them to calculated 

495 cloud top height on an individual grid cell basis. 

496 Assuming the density profile of the standard atmosphere, the total mass of 

497 lightning NO x was averaged over the EULINOX and STERAO storms. That amount of 

498 NO x distributed uniformly over a 200 km by 200 km region (typical of a global model 

499 grid cell) using the average midlatitude profile shown in Table 2 corresponds to a 

500 maximum increase in NO x mixing ratios of - 145 pptv between 7 and 9 km immediately 

501 following convection (~2 km higher than the maximum of the lightning NO x mass 

502 distribution). Because the lifetime of NO x increases with altitude, lightning NO x will be 
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503 converted to reservoir species such as PAN and HNO 3 more rapidly in the 7-9 km layer 

504 than at higher altitudes. As the time after convection increases, the maximum increase in 

505 NO x mixing ratios due to lightning would be seen at higher altitudes, which is consistent 

506 with the C-shaped profile of NO* typically observed in the troposphere. Downward 

507 transport from the stratosphere also contributes to the upper tropospheric maximum in 

508 observed NO x mixing ratios, while the maximum near the surface results from emissions 

509 from surface sources such as fossil fuel combustion and soil. 

5 1 0 The impact of assuming the vertical distributions of lightning NO x mass presented 

5 1 1 here has been calculated using NASA's Global Modeling Initiative (GMI) combined 

512 stratosphere-troposphere (Combo) CTM. The Combo CTM is detailed in Ziemke et al. 

513 [2006] and Duncan et al. [2007]. The tropospheric chemical mechanism includes 80 

5 1 4 species and 300 reactions to simulate 03 -NO x -hydrocarbon chemistry [Bey et al. 200 1 j . 

515 Lightning is assumed to produce 5 TgN yr' 1 . Horizontal distribution of lightning NO., is 

516 determined by the locations of parameterized deep convection (as indicated by the upper 

5 1 7 tropospheric values of cloud mass flux) in the meteorological fields from the Global 

5 1 8 Modeling and Assimilation Office's GEOS-4 Data Assimilation System that are used to 

519 drive the Combo CTM. Flash rates are scaled such that on a regional and monthly basis 

520 they match those from the OTD/L1S climatology (Allen et al., 2009). 

521 Two one-year simulations were produced by the Combo CTM. The first used the 

522 vertical profiles of lightning NO x from Pickering et al. [1998] while the second 

523 simulation used the vertical profiles of lightning NO x shown in Figure 12a and b for the 

524 midlatitude and subtropical regimes. Because no tropical thunderstorms were simulated 

525 with the cloud and chemistry models, a hypothetical tropical distribution (Figure 12c) 

526 was constructed by extrapolating the subtropical profile to a higher tropopause regime. 
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527 Test runs of the GMI Combo CTM revealed that this profile performed well in regions 

528 dominated by marine convection, but in tropical continental areas the Pickering et al. 

529 [1998] profile yielded better results for upper tropospheric ozone compared with 

530 ozonesonde data. Therefore, in tropical continental regions the Pickering et al. [1998] 

53 1 profile, modified by removing the boundary layer maximum and redistributing this mass 

532 between 4 and 1 1 km, was used. The tropical profiles are also provided in Table 2. 

533 Figure 13 shows the zonal mean change in NO x and O 3 in January and July when 

534 the profiles shown in Figure 12 are used instead of the Pickering et al. [1998] profiles. 

535 Because these profiles place less NO in the upper troposphere, NO x decreases above 1 0 

536 km with decreases over 1 00 pptv found over the southern hemisphere in January and over 

537 the northern hemisphere in July. The decrease in NO x in the upper troposphere, which 

538 results from using the modified vertical profiles of lightning NO x mass, also causes a 

539 small decrease in ozone throughout much of the troposphere. The largest decrease in 

540 ozone is ~10 ppbv in July at 15 km and 30°N. 

541 4.2. NO production 

542 A best-fit production scenario of Pic and Pcg has been estimated by comparing in- 

543 cloud aircraft observations with model output for the five storms where in-cloud aircraft 

544 observations were available. Figure 14 shows the production scenarios estimated for 

545 these five storms as well as the production scenario from Price et al. [1997], which was 

546 used in calculating the vertical profiles of LNO x mass presented in Pickering et al. [1998] 

547 and has been used in many global CTMs. In all cases, Pcg was estimated to be less than 

548 the 1100 moles per CG flash given in Price et al. [1997]. In addition, in all cases the ratio 

549 of Pic to P C g was greater than the commonly assumed value of 0. 1 presented by Price et 

550 al. [1997]. Over the five storms simulated, the average estimated Pcg was 500 moles NO 
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55 1 per flash or 7 kg N per flash (range of 360 to 700 moles per flash or ~5 to 9.8 kg N per 

552 flash). Assuming Pic/Pcg ratios in the middle of the estimated uncertainty ranges for the 

553 July 29 CRYSTAL-FACE, July 12 STERAO, and July 21 EULINOX storms yields an 

554 average Pic/Pcg ratio of 0.93, corresponding to 465 moles NO. The median peak current 

555 (16.5 kA for negative flashes and 19.8 kA for positive flashes, which account for 10.9% 

556 of the total) of the North American Lightning Detection Network (NALDN) presented in 

557 Orville et al. [2002] corresponds to a Pcg value of 508 moles NO when using the Price et 

558 al. [1997] relationship between peak current and energy dissipated, which agrees well 

559 with our estimate of 500 moles NO per CG flash. Therefore, the cases we have simulated 

560 appear to be representative of midlatitude and subtropical lightning. Assuming the 

561 average production scenario over the five storms presented, an average global IC to CG 

562 ratio of 3 (the same as estimated for the continental United States by Boccippio et al., 

563 2001) , and a global flash rate of 44 flashes s' 1 [Christian et al., 2003] yields a global 

564 lightning NO source of 8.6 Tg N yr' 1 . 

565 Estimates of the global lightning NO source have ranged from 2-20 Tg N yr' 1 

566 [IPCC, 2001], This range has been narrowed in recent years. IPCC [2007] recommends 

567 a 1 . 1 -6.4 Tg N yr' 1 range, though a more comprehensive review has specified the range to 

568 be 2-8 Tg N yr' 1 [Schumann and Huntrieser, 2007]. This review considered all of the 

569 methods of estimating LNO x production: global model estimates constrained by various 

570 observations, cloud-resolved modeling studies, aircraft measurements, laboratory 

571 experiments, and theoretical estimates. Our estimate of 8.6 Tg N yr' 1 lies near the upper 

572 end of the recently established range of global estimates. It is unlikely that our estimate 

573 is larger than values in this range due to the global extrapolation. Because of satellite 

574 observations of total lightning activity from instruments such as the Optical Transient 
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575 Detector (OTD; see Christian et al., 2003), uncertainty in the global flashrate has been 

576 greatly reduced. Perhaps more likely it is due to the fact that our average production 

577 scenario was calculated using data from only midlatitude continental and subtropical 

578 storms. No tropical thunderstorms were simulated in this analysis. Because 78% of 

579 lightning flashes occur between 30°S and 30°N [Christian et al., 2003], further 

580 investigation of the properties of tropical lightning flashes and their production of NO is 

581 needed. Based on analysis of data from the TROCCINOX experiment in Brazil, 

582 Huntrieser et al. [2008] suggest that tropical flashes (shorter, possibly due to low vertical 

583 wind shear) produce less LNO x than subtropical or midlatitude flashes (longer, possibly 

584 due to greater wind shear). If this hypothesis holds true throughout the tropics, inclusion 

585 of tropical events into our average production scenario will decrease this average val ue. 

586 Tropical event simulations are currently underway. 

587 Several modeling studies have assumed a production scenario similar to the mean 

588 scenario presented here where both IC and CG flashes produce on average approximately 

589 500 moles NO per flash. Hudman et al. [2007] simulated NO x over the United States 

590 during the 1CARTT campaign using the GEOS-Chem CTM and found that assuming a 

591 production per flash of 500 moles NO instead of the default GEOS-Chem value of 125 

592 moles NO improved the comparison with upper tropospheric aircraft observations. 

593 Cooper et al. [2006] simulated LNO x production and transport over North America 

594 during the same period using the FLEXPART model and assumed IC and CG flashes 

595 produce 460 moles NO (based on the DeCaria et al. [2005] results for the July 12 

596 STERAO storm) and a 2-day lifetime for upper tropospheric NO x . Resulting ozone 

597 enhancements were estimated with a box model. The results showed a good agreement 

598 between simulated NO x and DC-8 aircraft observations and between simulated ozone 
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599 enhancements and observations from commercial aircraft and ozonesondes. Jourdain et 

600 al. [2009] found that assuming a production per flash of 520 moles NO (rather than 260 

601 moles NO) in GEOS-Chem simulations of ozone over the United States substantially 

602 reduced the model’s bias relative to TES observations. Choi et al. [2008] compared 

603 results from a regional CTM with tropospheric NO 2 columns from OMI and tropospheric 

604 O 3 columns derived from a combination of OMI and MLS data and found that the model, 

605 assuming IC and CG flashes produce 500 moles NO, was able to reasonably reproduce 

606 the spatial distribution of the satellite observations. 

607 4.3 NO 2 in the vicinity of electrified storms 

608 NO 2 profiles appropriate for regions dominated by NO x from lightning are useful 

609 for improving satellite retrievals of NO 2 column amounts in such environments. The air 

6 1 0 mass factor used in converting satellite-measured slant columns to vertical columns of 

61 1 NO 2 is highly dependent on the NO 2 profile shape. NO 2 profiles from the version of the 

612 CSCTM that includes chemical production and loss have been used to examine the 

613 structure of this species in the CRYSTAL -FACE, STERAO, and EULINOX storms. 

614 Mean vertical profiles of NO 2 after convection are shown in Figure 15. The 

615 profiles were calculated over a 40 by 40 km area positioned in the convective core region 

616 of the model domain. The July 29 profile shows extremely elevated NO 2 mixing ratios 

617 exceeding 2 ppbv from 5 to 7.5 km, while the July 16 profile maximizes at only 0.7 ppbv 

618 at 9.5 km. The disparity in profiles is due to the large difference in flash rates in the two 

619 storms with over 4000 CG flashes recorded by the NLDN in the July 29 storm, and only 

620 301 CG flashes recorded in the July 16 storm. Of the midlatitude storms, the vertical 

621 profile of NO 2 from the EULINOX storm (with both high flash rate and pollution inflow) 
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622 is much larger than in either of the STERAO storms with a maximum of nearly 3 ppbv at 

623 5.5 km. 

624 In addition to NO 2 profiles, we have also used these simulations to estimate 

625 partial columns of NO 2 that a satellite may be able to observe in the upper portion of a 

626 convective cloud. Vasilkov et al. [2008] compared cloud-top pressures retrieved in the 

627 UV wavelengths by the Ozone Monitoring Instrument (OMI) aboard the Aura satellite 

628 with those retrieved by IR measurements from MODIS aboard Aqua. A comparison of 

629 both cloud-top pressure products with CloudSat radar observations indicate that OMI 

630 may be able to measure NO 2 as far down as 400 to 600 mb in the presence of deep 

63 1 convective clouds. Partial NO 2 columns were calculated for our storms of interest by 

632 averaging 2-km model output over 1 3 km by 24 km areas equivalent to OMI's nadir 

633 footprint and integrating the column from 400 mb to the tropopause. The results are 

634 shown in Figure 16. The July 29 CRYSTAL-FACE and July 21 EULINOX storms show 

635 the greatest degree of enhancement, which is consistent with the NO 2 profiles shown in 

636 Figure 15. In the July 29 CRYSTAL-FACE storm, the region of highly elevated NO;> 

637 columns extends over 1 00 km. Table 3 shows the sensitivity of the peak and background 

638 NO 2 columns to the choice of column depth. The peak values indicate the maximum 

639 partial NO 2 column amounts over an OMI grid box area (13 km by 24 km) in the storm 

640 anvil or core regions, while the background values indicate the mean partial NO 2 column 

641 amounts outside of the storm region. These data should be useful in providing a range of 

642 values of tropospheric column NO 2 that may be compared to retrievals over convective 

643 clouds from such instruments as OMI. 

644 5. Conclusions 
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645 Simulations of six midlatitude and subtropical thunderstorms occurring during 

646 four field projects have been conducted using the CSCTM, which includes a 

647 parameterized source of LNO x . In order to estimate LNO x production per flash in each 

648 storm, different scenarios of Pie and Pcg were specified in the model, and the results 

649 compared with in-cloud aircraft observations of NO x . By comparing column mass and 

650 probability distribution functions of the observed and simulated storms, the most 

65 1 appropriate production scenario was estimated for each storm. The results suggest that 

652 Pic may be nearly the same as Pcg- A frequent assumption that Pic is equal to one tenth 

653 Pcg resulted in a significant underestimation of LNO x in all five simulations where anvil 

654 observations were available. We echo the recommendation of Ridley et al. [2005] for use 

655 of comparable values of Pcg and Pic. The Ridley et al. [2005] recommendation was 

656 partly based on one of the cases presented here. Our similar results for four additional 

657 cases strengthen this recommendation. The mean Pcg obtained from the five case 

658 simulations is 500 moles NO per flash (range 360 to 700), which, when extrapolated 

659 globally, yields an estimate at the high end of the currently accepted range. This may 

660 result from the absence of tropical flashes in this analysis. 

661 Vertical profiles of the percentage of LNO x mass in each 1-km layer after the 

662 convection show very little LNO x mass near the surface with the majority of LNO x 

663 remaining in the mid- and upper troposphere in a “backward C-shaped” profile. Global 

664 and regional CTMs that have adopted C-shaped vertical profiles of LNO x mass may be 

665 underestimating the amount of LNO x in the mid- troposphere and overestimating the 

666 amount near the surface. Global models which are changed to represent equivalent pet- 

667 flash production by IC and CG lightning may require a change of the global LNO x source 

668 strength to reasonably reproduce NO x observations in the middle and upper troposphere. 
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669 Changes in the vertical placement of LNO x in CTMs may significantly alter distributions 

670 of species such as O 3 and OH. 

671 
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937 Table 1 . Calculated column mass of N in NO x in the July 29 CRYSTAL-FACE storm. 

938 assuming an IC/CG ratio of 5 

939 


Pcg (moles NO 
per flash) 

Pic/PCG 
(IC/CG = 5) 

Pic/PcG 
(IC/CG = 2) 

Column Mass 
(gNm- 2 ) 
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6.9 x 10' 4 

590 

0.75 

1.9 

8.3 x 10' 4 

590 

1.0 

2.5 

1.0 x 10' J 
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940 Table 2. Average Profiles of LNO x Mass in Percent 

941 


Altitude 

Range, 

km 

Subtropical 

Tropical 

Midlatitude Continental 

Tropical 
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945 Table 3. Partial NO 2 columns (10 14 molecules cm' 2 ) from the tropopause to 400, 500 am 

946 600 mb 

947 

948 Tropopause to 400 mb Tropopause to 500 mb Tropopause to 600 mb 

949 Storm Peak Background Peak Background Peak Background 
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Figure 2. Pdfs of simulated and observed NO x at 12.5 km (a) and 13 km (b) for the July 
29 CRYSTAL-FACE storm. Solid (dashed) line shows observed values assuming clear 
sky (cloud enhanced) photolysis rates. Dashed (dot-dashed) line shows simulated values 
assuming Pic/Pcg of 0.5 (0.6). 
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984 


Relationship between P CG and P iC /PcG 



986 Figure 3. Relationship between Pcg and the Pic/Pcg ratio necessary to match the column 

987 mass of N computed from aircraft observations. Solid (dashed) line shows the 

988 relationship when an IC/CG ratio of 2 (5) is assumed. Plus signs indicate the production 

989 scenarios when Pcg is assumed to be 590 moles NO per flash based on observed mean 

990 peak current. 
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99 1 (MBS: COU/IOE3 

992 

993 Figure 4. Vertical cross-section of simulated N0 X in the July 29 CRYSTAL-FACE storm 

994 (assuming P C q = 590 moles NO, Pic = 354 moles NO, and an IC/CG ratio of 5) at 240 

995 minutes. 
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%/km of Total Lightning N Mass 

Figure 5. Vertical distributions of the percentage of LNO x mass per kilometer injected 
into the cloud and following convection for the July 29 CRYSTAL-FACE storm 
(assuming Pcg = 590 moles NO, Pic = 354 moles NO, and an IC/CG ratio of 5). 
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Figure 7. Vertical cross-section of simulated NO x through the core of the July 16 
CRYSTAL-FACE storm assuming Pcg = 700 moles NO and Pic = 630 moles NO at 180 
minutes. 


44 




0 5 10 15 20 25 

%/km of Total Lightning N Mass 


1014 

1015 

1016 Figure 8. Vertical distributions of percentage of LNO x mass per kilometer following 

1017 convection for two simulated subtropical storms. 
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Figure 9. Vertical distributions of the percentage of LNO x mass per kilometer following 
convection for (a) the July 16 and (b) the July 29 CRYSTAL-FACE simulations 
assuming IC to CG ratios of 2 and 5. 


46 





1026 

1027 

1028 

1029 

1030 

1031 

1032 


NO x (ppbv) 180 minutes y = 68 



CD 

~o 

Z3 

H-J 

- P 

<c 


Distance (km) 


km 



Figure 1 0. Vertical cross-section of simulated NO x through the core of the southern ceil 
of the July 10 STERAO storm assuming Pcg=390 moles NO and Pic=234 moles NO at 
180 minutes. 
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1032 

1033 

1034 Figure 11. Vertical distributions of percentage of LNO x mass per kilometer following 

1035 convection for four simulated midlatitude continental storms. 
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Figure 12. Average vertical distribution of percentage of LNO x mass per kilometer 
following convection (solid) for the (a) subtropical and (b) midlatitude continental 
regimes. Dashed line shows polynomial fit. Midlatitude continental profile from 
Pickering et al. [1998] (dash-dot) is also shown in (b). Hypothetical tropical marine 
profile (c) is based on extrapolating the subtropical average profile to a higher tropopause 
regime while the tropical continental profile was constructed using the Pickering et al. 
[1998] profile with the boundary layer maximum removed and that mass redistributed 
into layers from 4 to 1 1 km. 
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1053 Figure 13. Zonal mean change inNO x (pptv) and O 3 (ppbv) in January and July resulting 

1054 from using the lightning NO x profiles in Figure 1 1 instead of the profiles from Pickering 

1055 etal. [1998], 
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1056 

1057 

1058 Figure 14. Estimated lightning NO x production scenarios for the July 16 CRYSTAL- 

1059 FACE (C, 7/1 6), July 29 CRYSTAL-FACE (C, 7/29), July 10 STERAO (S, 7/10), Julv 

1060 12 STERAO (S, 7/12), and July 2 1 EULINOX (E, 7/2 1 ) storms. The bars in the 7/29 

1061 CRYSTAL-FACE, 7/12 STERAO, and 7/21 EULINOX storms indicate the estimate 

1 062 uncertainty. The Price et al., [1997] and Fehr et al. [2004] production scenarios are 

1063 indicated by asterisks, as is the estimated value of P C g calculated assuming the NALDN 

1 064 median peak current from Orville et al. [2002]. 
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(a) Subtropical N0 2 




Figure 15. Average vertical profiles of NO 2 mixing ratio in the convective core region 
following simulated subtropical storms (a) and midlatitude storms (b). 
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1072 

1073 Figure 16. Partial NO 2 columns (10 14 molecules cm' 2 ) from the tropopause to 400 hPa 

1 074 for CRYSTAL-FACE, STERAO, and EULINOX simulated thunderstorms. 

1075 


53 


